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ABSTRACT 

Sulfur abundances derived from optical emission line measurements and ion- 
ization correction factors in planetary nebulae are systematically lower than ex- 
pected for the objects' metallicities. We have carefully considered a large range 
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of explanations for this "sulfur anomaly", including: (1) correlations between the 
size of the sulfur deficit and numerous nebular and central star properties; (2) 
ionization correction factors which under-correct for unobserved ions; (3) effects 
of dielectronic recombination on the sulfur ionization balance; (4) sequestering 
of S into dust and/or molecules; and (5) excessive destruction of S or production 
of O by AGB stars. It appears that all but the second scenario can be ruled out. 
However, we find evidence that the sulfur deficit is generally reduced but not 
eliminated when S"*"^ abundances determined directly from IR measurements are 
used in place of the customary sulfur ionization correction factor. We tentatively 
conclude that the sulfur anomaly is caused by the inability of commonly used 
ICFs to properly correct for populations of ionization stages higher than S"^^. 

Subject headings: ISM: abundances, abundances, planetary nebulae: general, 
stars: evolution 



Introduction 



The sulfur anomaly was first identified by iHenry et al. I (120041 hereafter HKB04) and 
refers to the systematically lower sulfur abundances exhibited by Galactic planetary nebulae 
(PNe) compared to H II regions of the same metallicity, where metallicity is gauged by 
oxygen abundance. Specifically, HKB04 computed S abundances by observing S+ and S"*"^ 
abundances directly and then applying a sulfur ionization correction factor (ICF) to correct 
for the contribution of unobserved sulfur ions to the total abundance. HKB04 were unable 
to establish an exact cause but speculated that the problem is likely brought about by the 
use of a sulfur ICF which fails to correct adequately for unobserved ions. 



Th e sulfur anomaly has also been reco g nized and disc u ssed more recently bylBeriiard-Salas 



f2006h. 


Pottasch & Bernard-Salas 


( 


2006). 


Shaw et al. 


Green et al. 


(2011 


), and 


Rodriguez & Deleado-Inelada 


on its cause has developed. 


Pottasch & Bernard-Salas 



mm 
(Eoil. 



Kwitter fc Henrv I ( 1201 If ) 



but still no clear consensus 



mm and iBernard-Salas I mm 



used ISO measurements in the IR to obtain S^^ ion abundances of 26 PNe. They then 
obtained the total sulfur abundance by computing the sum of the S+, S"*"^, and S"*"^ abun- 
dances for each object, using abundance values found in the literature for the first two ions. 
Their resulting S abundances still fell below expected levels, and they suggested that the 
missing gas-phase sulfur might b e sequ estered into dust. (We address this idea below in Sec- 



tion HI) Similarly, I Shaw et al. I ( 120101 ) used both deep ground-based optical spectra along 
with Spitzer IR me asurements to determine su l fur ab unda nces in 14 PNe i n the Small Mag- 
ellanic Cloud. Like jPottasch fc Bernard-Salas I (|2006[ ) and iBernard-Salas I (j2006[ ) they found 
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that the direct observation of S"*"^ did not solve the problem. They also claimed that the sul- 
fur abundance problem worsened for their three most highly ionized objects. This appears to 
be direct evidence suggesting that ionizatio n stages above S"*"^ are significan t ly po pulated yet 
unaccounted for by a standard sulfur ICF. [Rodriguez fc Delgado-Inglada I (1201 if ) computed 
models of PNe with various physical conditions and then subjected the output spectral line 
strengths to the same analysis as real observed measurements of PNe. They found that the 
S abundance inferred from the predicted line strengths can exceed the input abundance by 
0.05 dex or fall below it by 0.3 dex and s tressed that c a ution shou ld be used when using an 
ICF to determine S abundances. Finally, iGreen et al. I ( 120111 ) and iKwitter fc Henry I (120111 ) 
discussed the sulfur anomaly but offer no solutions to the problem. 



The sulfur anomaly can be seen clearly in Figured], where we have plotted e(^SjLJ versus 
e(0) for a large sample of PNe in the disks of the Milky Way and M31 galaxies (open 
symbols) as well as a combined sample of blue compact galaxies (BCG) and Galactic and 
extragalactic H II regions (hereafter referred to as H2BCG; filled circles) for comparison 
purposes. Fig. [2] is an analogous plot of e(Ne) versu s e(0). In both plots the a bundances for 



the M ilky Way PN sample were taken from HKB04. lMilingo et al. I (120101 ). and lHenry et al. 



(120 lOl ). where data for all objects were obtained, reduced, and measured by the same team 
of astronomers. In addition, all elemental abund ances were recently recom puted using an 
updated version of the ELSA abundance package ( Johnson fc Levitt! bood ) and the results 
compiled into a database by Dr. Karen Kwitter and her students at Williams College. Thus, 
the final abunda nces are considered to be homogeneous. Likewise, the M31 PN abundances 

J2012h . 



were taken from Kwitter et a. 



PNe from lLeisy fc Dennefeld 



In Fig. [2] we have added measurements of 70 LMC 
( 20061 ). Abundances for the H2BCG comparison samp l e wer e 



taken directly from the literature. The BCG data were tak en fromllzotov fc Thuan I (119991). 



For t 



le H I I regions we u s ed ab undances in MlOl from iKennicutt. Breso^ 



in. fc Garnett 



(2003) , from IShaver et al. I (jl983l ) for three objects in the MWG, and from lEsteban et al. 

( 2004 ) for the Orion N ebula. The solid line is a least squares fit to the H2BCG samplr^ 
Solar abundances from lAsplund et al. I (120091 ) are shown with dashed lines. 



The difference in behavior between the PN sample and the sample of H2BCG in Fig. [T] 
is striking. For the latter group, the objects clearly fall along a narrow linear track which 



ie(X) = 12 + log(X/H) 

^To check on the accuracy of the H2BCG abundances, we compiled the emission hne strengths employed 
by the various authors and recomputed the abundances using ELSA. Wc found the ELSA oxygen abundances 
were generally 1-2% higher than the published values, while the ELSA neon and sulfur abundances were 
each roughly 10% higher than the published numbers. Because of the small differences, we proceeded by 
using the original abundance values found in the literature. 
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extends over roughly two orders of magnitude in metallicity. This lockstep behavior between 
S and O is expected, since these two elements are alpha elements, synthesized under similar 
conditions in massive stars. Assuming the existence of a universal initial mass function, 
the lockstep behavior should be galaxy- independent. As H II region and BCG abundances 
represent element levels in the ISM of their host systems, their behavior exhibited in Fig. [1] 
is not surprising. 

At the same time, PN abundances reflect not only the level of an element present in the 
progenitor star at the time that it formed but also any changes to those abundances brought 
about by nuclear processes occurring within the star during its lifetime. In the specific cases 
of S and O, there is currently no strong evidence, except perhaps at very low metallicity, 
that PN progenitor stars alter their initial levels of these two elements (see the extended 
discussion on this topic in section [S] below), and so their abundances directly reflect the 
amounts of S and O present in the local ISM when the star formed. Under this assumption, 
we would expect the lockstep behavior seen in the case of the H2BCG objects to be the same 
for PNe. As we see in Fig. [H it is not. 

Interestingly, the offset in S abundances between PNe and H2BCGs is not duplicated 
for neon abundances shown in Fig. [2l i.e., we see no evidence for an offset between the two 
object types, although PNe continue to show more scatter. Therefore the differences in PN 
behavior demonstrated in Figs. [1] and [2] suggest that the offset in Fig. [1] is related to sulfur 
and not oxygen. Analogous plots of S versus either Ne or Ar show that S is always offset from 
the H2BCGs by roughly the same amount and in the same direction, further emphasizing 
that the problem resides with sulf ur. A similar conclus ion is drawn when one inspects the 



plots of Ne, Ar, and S versus O in iMilingo et al. I ( l2010l ). 



In what follows, we carefully examine a large number of potential explanations for the 
sulfur anomaly. In section [2] we test for correlations between the magnitudes of the vertical 
offsets of PN sulfur abundances from the H2BCG track in Fig. [1], i.e., sulfur deficits, and 
numerous independent parameters related to the nebula and/or central star. In section [3] 
we check on the robustness of the sulfur ionization correction factor (IGF) using standard 
photoionization models, tests of the relevance of dielectronic recombination, as well as direct 
measurements of the S+^ abundance available in the literature. Section H] looks at the 
possibility that S is sequestered in either dust or molecules, while section [5] explores the 
possibility that the sulfur anomaly results from stellar nucleosynthesis in AGB stars. Our 
summary and conclusions are provided in section [6l 
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2. Correlation Tests Involving Nebular and Stellar Properties 

Table [1] lists the results of correlation tests performed on our PN database. We define 
the sulfur deficit SD = e(Sexp) — e(Sicf), where e(Sexp) is the expected S abundance, i.e. the 
abundance that a PN of a given O abundance would possess if its S abundance placed it on 
the H2BCG track, while e(Sicf) is the S abundance of a PN inferred with the use of an ICF. 
Thus the more that a PN is displaced below the H2BCG track the more positive the SD 
value. For each test the sulfur deficit was plotted against one of the 24 separate independent 
parameters listed in column 1 for the number of objects in column 2. The variables tested, 
from top to bottom in column 1, are: reddening factor, log of the ion ratio, oxygen ionization 
correction factor, sulfur ionization correction factor, log of the nebular diameter, central star 
mass, log of the central star's luminosity, neon deficit (defined in the same way as the sulfur 
deficit), log oxygen abundance, log of the ion ratios (4 rows), log of the difference between 
a PN's measured O abundance and the ambient interstellar abundance at its location along 
the disk, the PN's galactocentric distance along the disk, electron temperatures derived from 
the lines indicated (next five variables), electron densities based on the lines indicated (next 
two variables), and excitation class based up on the indicated emission l i nes an d employing 



expressions for excitation class developed by iDopita &: Meatheringham I (119901 ). 



Column 3 displays the value of the Pearson correlation coefficient, while column 4 lists 
the correlation probability, i.e., the probability that the correlation could have resulted by 
choosing the same number of objects from a completely uncorrelated parent population|j. 
Defining r=0.6 as the minimum value for a likely correlation (our justification of this is ex- 
plained in Appendix A), all but one of the tests indicate the absence of a correlation. The 
exception is the relation between the SD and log(S+'^/S"''^). The increase of this ratio with 
the SD suggests that the deficit is due in part to the failure of sulfur ICFs to correctly ac- 
count for the higher ionization stages of sulfur. Interestingly, we see no significant correlation 
between the SD and either the Ar"'"^/Ar"'~^ or Ar''''^/Ar+^ ratios. The ionization potentials of 
S+2 and S+^ are 34.8 eV and 47.2 eV, while those of Ar+2 and Ar+3 are 40.7 eV, and 59.8 eV, 
respectively. Since we see no significant correlation for Ar+^/Ar+^, this suggests that there 
should be little sulfur in stages above S"*"^. If the sulfur anomaly is the result of underes- 
timating the population of ionization stages above S"*"^, it appears that the only significant 
stage in this case would be S+'^. We conclude that, with the exception of log(S"'"^/S+^), none 
of the other factors listed in Table [T] is significantly related to the sulfur anomaly. We should 
point out, however, that the question of a link between sulfur deficit and PN morphology 
was not explored. Such a study should perhaps be undertaken in the future. 



^We employed the pearsn subroutine found in Press et al. ( 20031) 



- 6 - 



3. Evaluating the Effectiveness of the Sulfur Ionization Correction Factor 



The ionization correction factor is the ratio of the elemental abundance to the total of 
the observed ionic abundances. It corrects for the unobservable ions of an element. In the 
case of sulfur, normally only the S"*" and S"*"^ ions are directly observable, although we know 
from photoionization models that higher stages, e.g., S"*"^, are present in PNe in significant 
quantities. ICFs are usually based upon nebular models from which the abundances of the 
unseen ions can be estimated. It is therefore natural to speculate that the observed PN 
sulfur deficits are due to a formulated ICF for sulfur which incorrectly adjusts for the unseen 



ions. 



Various formulas for the sulfur ICF have been deve 



oped over the l ast four decades, in- 



cludin g thos e of iPeimbert fc Costero I(ll969l ). IStasiriska I (119781 ) , [French I (Il98ll ) , iKingsburgh fc Barlow 
(I1994J ). and iKwitter fc Henry I ( 120011 ). With the advent of ISO and Spitzer, however, the 
S~^^ ion has been observed directly in oyer 30 PNe in which oxygen abundances were also 
reported by iPottasch fc Bernard-Salas I (120101 ). These observations have lessened the need 
of an ICF in these PNe and allowed for the more direct computation of the S abundance, 
although the sulfur deficit problem has not been discussed extensively in the context of these 
observations. 

In this section we check to see if photoionization models, from which ICFs are derived, 
reproduce the sulfur deficits observed in three sample PNe. We also test the relevance of 
dielectronic recombination and the possibility that incorrect rates could produce flawed ICFs. 
Finally, we report on the implications that direct IR measurements of S~^^ have for the sulfur 
deficits. 



3.1. Photoionization Models 



3.1.1. Standard Models 



The abundance analysis by HKB04 of their PN sample included the calculation of 
photoionization models of three PNe representing a large range in the sulfur d eficit, i.e.. 



IC 4593, Hu 2-1, and NGC 3242. The program CLOUDY (IFerland et al. Ill998l ) was used 



along with input blackbody stellar SEDs to model each object's measured line strengths. In 
the process, elemental abundances, which comprised part of the model input information, 
were adjusted along with other parameters such as gas density and central star temperature 
and luminosity in order to produce a list of predicted line strengths which closely matched 
the observed values. [Since then, the sulfur abundance of NGC 3242 has been revised and 
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a new model of that PN was recently calc ulated b y Mat thew Hosek at Williams College 



using a central star SED of 91,000 K from iRauch I (Il997l ). Results of the new model are 



used in this discussion.] Table |5] provides a comparison of empirical and model O and S 
abundances along with associated sulfur deficits. We see that the model abundances for O 
and S along with the resulting sulfur deficit agree well with the analogous values inferred 
from observations for each of the three PNe used in this analysis. 

Clearly the fact that models and observations agree closely suggests that either standard 
models should not be used to derive sulfur ICFs, or the sulfur ICF is correct and there is 
a real gas-phase sulfur abundance deficit in these nebulae. In the former case, some S may 
be sequestered in ionization stages whose populations are underestimated by the models for 
reasons that are unknown at this time. Next, we consider possible effects of dielectronic 
recombination on the ionization structure of S in model nebulae, and by extension, the ICF 
for sulfur. 



3.1.2. Dielectronic Recombination 



Many of the sources of atomic data used in the photoionization code CLOUDY are 
summarized in lFerland et al. I f ll998l ) with updates on the Cloudy web site www.nublado.org. 
Could uncertainties in the atomic data produce the sulfur anomaly discussed here? 

The ionization balance of hydrogen and helium are quite reliable, since the ionization and 
recombination of these elements is dominated by valence shell photoionization and ra d iative 



recombination. These processes are discussed, for instance, in lOsterbrock fc Ferland I (12006 
Chapter 2). 

The situation is complex and more uncertain for the heavy elements. Ionization for a 
PN is predominantly by valence shel l photoionization, although i nner shell photoionization 
and Auger decay can be significant (lOsterbrock fc Ferland I l2006l . Chapter 11). The heavy 
elements primarily recombine through dielectronic recombination, a process where a free 
electron is captured into autoionizing levels. These are levels wit h two excited electrons 



and with energies above the ionization limit of the recombined atom ( lOsterbrock &: Ferland 



2006 



Appendix A4). These rates of dielectronic recombination (DR) are sensitive to the 
energies of the autoionizing levels because they determine whether free electrons can reach 
those levels. In practice, accurate DR rates require experimental autoionization level ener- 
gies. 



For high gas temperatures, we use the recombination data set collected by Nigel Badnell 
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for both radiative and dielectric recombinatioi]0. DR rates are available for elements of the 
second row (O and Ne are important for optical spectra of PN) but only for relatively high 
ions of the third row. In particular, experimental energies are not available for the low ions 
of sulfur studied in this paper, and their sulfur DR rates are not given by Badnell. Instead, 
for these ions we use estimates of the rates based on the simple averaging approach described 



by lAli et al. I (jl99ll ). That is, the Badnell rates are joined on to the mean of charge-specific 



CNO rates at lower temperatures. 

It is impossible to say how uncertain these estimates might be. Comparing with other 
ions this method generally works to within a factor of two. Could such uncertainties affect the 
ionization balance and line emission enough to cause models to misrepresent what happens 
in real nebulae? We vary the estimated DR rate by these amounts to find out. 

In the present case, we investigated the effects of varying the DR rate on sulfur emission 
line strengths as a way of identifying a mechanism which could profoundly affect the ICF of 
sulfur in a way that would explain the observed sulfur deficits. For example, overestimated 
rates would affect the sulfur ionization balance by shifting sulfur into lower stages, causing 
an ICF based on these models to underestimate the sulfur elemental abundance. An opposite 
shift would result if rates used are below their actual values. 

We therefore computed three grids of photoionization models using CLOUDY 08.01, 
where input parameter ranges for each grid were designed so that the output predicted line 
strengths spanned the space of observable line ratios for the PN sample. Three model input 
parameters were varied: the SED of the central star, the overall metallicity, and the nebular 
ga s density. For s tellar temperatures between 30,000 and 50,000 K we used the Atl as SEDs 



bv iKurucz I ( Il991l ). while for temperatures above 50,000 K we employed the SEDs of iRauch 
( 1l997l ). For each grid, input stellar temperature ranged from 30,000-158,000 K, total gas 
density from 10^-10^, and metallicity fro m 0.25-2.5 time s solar. The DR rates in the first 



grid were the same as those prescribed by lBadnell I (Il99ll ). For the other two grids, the rates 



were halved or doubled. All models were radiation bounded. 

Results for the first grid are shown in Fig. |21 where for each model we computed the pre- 
dicted line strength ratio of ([S III] A9523+A9069)/[S II] A6724 and ([O III] A5007+A4959)/[O II]A3727. 
Predicted ratios for models employing the Atlas or Ranch SEDs are shown with open squares 
and circles, respectively, while the line ratios for the same MW disk PN sample shown in 
Fig. [1] are shown here with filled circles. First, we note that these models mostly span 
the space occupied by the observations. The exceptions appear to be the objects with low 



''These rates are available at the following websites: |http:/ /amdpp. phys.strath.ac.uk/tamoc/DR/ (dielec- 
tronic recombination) and http://amdpp.phys.strath.ac.uk/tamoc/RR/ (radiative recombination). 
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[S ni]/[S II] ratios along with a few with high [O III]/[0 II] ratios. In the first case, the 
observed reduced sulfur ratios are likely the result of partial telluric absorption of the near 
IR lines at AA9069,9532. The second group is more difficult to characterize. However, in 
comparing a number of their observed properties with those of the majority of the sample, 
all members of the high [O III]/[0 II] group are likely to be matter bounded. For example, 
compared with the rest of the sample, these objects all possess relatively low observed values 
of [O I] A6300, [O II] A3727, and [N II] A6584, i.e., lines typically reduced in strength relative 
to those produced by higher ionization stages when the nebula is matter bounded. 

Figure m demonstrates the effects of multiplying the Badnell DR rates by factors of 0.5 
and 2. Here we plot results for the same line ratios as in Fig. [3] for all three grids. Only the 
models employing the Ranch AGB SEDs are shown. The figure clearly shows the effects of 
changing the rates. Typically the model points shift up by about 0.1 dex when the rate is 
halved and down by a similar amount when the rate is doubled, the directions of the shifts 
being in the sense that is expected. The thermostat effect, where changes in the abundance 
of an ion do not linearly affect its emission lines due to energy balance, accounts for the 
relatively modest change in line emission, despite large changes in the assumed DR rates. 
Shifts are similar for those models based upon the Atlas SEDs. These changes in DR rates, 
representing reasonable uncertainties in their actual values, correspond to shifts in the sulfur 
ICF of about 13% or 0.05 dex above or below the value obtained before DR scaling. Clearly, 
these small changes in the ICF cannot explain the sulfur anomaly. 



3.2. Direct Infrared Observations of S"*"^ 



All of the sulfur abundances shown in Fig. [T]result from applying an ICF to the observed 
abundances of S+ and S+^. Does replacing the ICF by directly measuring S~^^ reduce or 
eliminate sulfur deficits? 

Sulfur abundances of 32 Galactic PNe determined thro ugh direct observations of S''"'^ 
have been summarized by iPottasch fc Bernard-Salas I (120101 ). Their study made use of the 
ISO and Spitzer telescopes to measure t he strength o f the J S IV] 10.5/im line. Similar 



I 



me asurements using Sp itzer were made by IShaw et al. I ( 120101 ) of six PNe in the SMC and 



by lHenry et al. I (120081 ) of DdDm-1, located in the Galactic halo. 



In Table |3] we provide S abundance results which have been derived using both the ICF 
method (Sicf) and the IR method (Sir). For each object identified in the first column we list 
in columns 2-5 the values for e(0); e(Sexp), i-e., the e(S) value on the H2BCG track for the 
given oxygen abundance; e(Sicf); and SDicf. For the last quantity we have added the subscript 
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to indicate that it is based on the standard ICF abundance method. Columns 6 and 7 give 
e(Sir) resulting from the IR measurements and the associated sulfur deficit SDir based upon 
values in column 3. Column 8 provides a measure of the difference in S abundances derived 
from the two methods. Finally, values in column 9 represent the amount of an unidentified 
ion of S such as S° or S"*"^, relative to Sexp, needed to explain the difference between Sir and 
Sexp, i.e., SyScxp = 1 - log~^(-SDh.). 

Fig. O shows a visual comparison of the e(Sir) and e(Sicf) values in Tabled where Milky 
Way disk objects appear as filled circles, SMC objects as filled diamonds. The dashed line 
indicates the one-to-one correspondence. With a few exceptions there is a clear tendency (21 
of 29 objects) for the IR abundances to exceed the corresponding ICF abundances. Of those 
21 PNe the median positive vertical distance from the line is 0.17 dex or a factor of 1.5 times 
the Sicf value. This suggests that when S"'"'^ observations replace the ICF method, half again 
as much S, on average, is inferred. This result is qualitatively consistent with the positive 
correlation we found in § 2 between the SD and log(S"'~^/S~'"^). For those eight objects falling 
below the line, apparently the ICF overcorrects for unobserved sulfur ions, as would be the 
case if these PNe were matter-bounded. 

Next, if the sulfur deficit is caused by the failure of the ICF method to properly correct 
for S+^, then there should be a direct correlation between the sulfur deficit SDicf and the 
quantity e(Sir) — e(Sicf). Values listed in columns 5 and 8 of Table [3] are plotted against each 
other in Fig. [HI For this analysis we focus on those objects with positive values for SDjcf , i.e., 
those to the right of the vertical dashed line, ignoring the one outlier near the vertical line at 
the bottom of the plot. Here we see by inspection some evidence for a positive correlation. A 
regression analysis of the 21 points provides the least squares fit shown by the solid line. The 
slope of the line is +0.34±.13 with an intercept of -0.01±.06 dex. However, the correlation 
coefficient r = +0.52, with a probability of less than 1% that the same set of points could 
have come from an uncorrelated sample. As explained in Appendix A, most statisticians 
require a value of r > 0.6 for a true correlation, and so we conclude that the apparent trend 
in Fig. [6] is not robust enough to confirm a relation between values in columns 5 and 8. 

Approaching this line of investigation reg arding S"*"^ from another angle, ho wever, we 



have added the S and O abundance data from iPottasch fc Bernard-Salas I (120101 open dia- 



monds) and lShaw et al. I (120101 . open squares) to the original plot in Fig. [T]to produce Fig. [71 
Compared with the original PN sample of open circles, we now see a reduced tendency for 
these objects to fall systematically below the H2BCG track. Also, the vertical scatter is 
much closer to what we see in the case of neon in Fig. [2] and perhaps can be attributed sim- 
ply to natural abundance scatter. Thus, the direct use of the S"*"^ abundances to determine 
total S abundances seems to solve at least part of the problem related to the sulfur anomaly. 
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as originally anticip a ted b y HKB04. This finding is also consistent with the model results 



of iGoncalves et al. I ( 120111 ). who use the 3-D photoionization code Mocassin to show that 
as the temperature of the PN cen tral star increases, the sulfur i on population shifts upward 
and the ratio of the sulfur IGF of iKingsburgh fc Barlow I fll994l ) to the 'true' value inferred 



from the ionization populations in the model output declines. 

Shifting attention now to the points in Fig. |6] with negative SDicf values, we see that 
there is no apparent relation between the parameters plotted. If we speculate that the sulfur 
anomaly is the result of an unknown nebular property that populates unseen ionization stages 
of S beyond what we find in the models discussed above, then we can conclude that this 
same property is not relevant for these seven objects. For them, other factors are influencing 
the values of Sir — Sicf, but they appear to be unrelated to the sulfur anomaly. 

We turn now to the interpretation of the values posted in the last column in Table|31 i.e., 
SyScxp [= 1 ~ log~^(~SDir)]. For a given oxygen abundance, if we take Sexp to be the sulfur 
abundance of a PN falling on the H2BCG track, then the value in column 9 is the fraction of 
that total abundance that represents the difference between it and the sum of S"*", S"*"^, and 
g+3 measurements, i.e.. Sir. (The six objects with negative SDir values are irrelevant here.) 
The average value of SyScxp is 0.43±.24. This is rougly the level of the estimated natural 
scatter inferred from neon abundances of about 0.2 dex or 60% (see Fig. |2]). Generally, then, 
replacing the IGF with direct S"*"^ measurements appears to reduce sulfur deficits in most 
PNe to the level of expected scatter in the S abundance. However, we note that these objects 
still fall systematically below the H2BGG track, and thus the deficits of PNe are not entirely 
explained by empl oying the direct measu reme nt of S^^. This conclusio n is consistent with 



results reported by iBernard-Salas I (120061 ) and iPottasch fc Bernard-Salas I ( l2006l ). 



Values in column 9 surpassing this scatter level would suggest that additional ions of 
S such as S° and/or S"*"^, are present in significant amounts in these objects. We note that 
an unexpe c tedly high abundance of ions such as S+^ in PNe is predicted by nebular models 
by iJacob I (120111 ) in which the density profile of the nebula is non-uniform. Our standard 
GLOUDY models suggest that S"*""^ has a permitted line S V A1198 corresponding to the 3s3p 
^P to 3s^ ^Sq transition. If we assume an ion fraction of 10% for S"^^, then, extrapolating 
from these same model predictions, the S V line will be well below detection levels. Note, 
however, that considerations of ion ratios and ionization potentials in § 2 appear to conflict 
with the idea of significant abundances of ions above S"*"^. 



-VI- 



Dust and Molecules 



In studies of the ISM, deficiencies in the abundances of element s in atomic (or i onized) 



gas i s usually assumed to be due to depletion into dust grains (e.g. IWhittet Ill992l : iKrugel 



20081 ). Consequen tly, the observed sulfur defic it in PNe has been attributed to sequestration 



into dust grains ( Bernard-Salas et al. II2008I ). The chemical composition of dust in PNe 



should reflect that of their parent Asymptotic Giant Branch (AGE) stars. The atmospheres 
of AGB stars are expected to change as these stars evolve, due to convective dredge up of 
carbon produced in the He-burning shell. The amount of carbon relative to oxygen (the 
C/0 ratio) is critical in determining which types of dust and molecules are present around 
an AGB star. The formation of extremely stable CO molecules will consume whichever of 
the two elements is less abundant, leaving only the more abundant element available for 
dust formation. Stars start their lives with the cosmic C/0 ratio of f5i0.4 and are therefore 
oxygen- rich. In about a third of AGB stars, enough carbon will be dredged up to make C/0 
> 1 and therefore carbon will dominate the chemistry around these stars. 

Sulfur-bearing solids (e.g. MgS, CaS and FeS) are pr edicted to form in carbon-rich 



environments (e.g. iNuth et al.l 119851 : iLodders &: Fegleyl 119951 ) and have been tentatively ob- 
served in carbon-rich AG B stars, PPNe and PNe, e.g. the "30 Atm " feature oftentimes 



associated with MgS (e.g. 


Goebel & Moselev 


1985: 


Nuth et al. 


Beeemann et al. 


1994; 


Honv. Waters. &: Tielens 


2002). 



In oxygen-rich environments sulfur-rich grains are precluded by the abundant free oxy- 
gen, which will tend to oxidize any sulfides that form, releasing the sulfur back to the gas 
phase. Therefore, we should expect sulfur depletion in the gas phase to be highest in carbon- 
rich nebulae. Fig. [S] shows that the distribution of sulfur deficits does not correlate with the 
the carbon abundance (the determination coefficient r^ is 0.202 or r=0.45; see Appendix A 
for a brief discussion of correlation coefficients.) The figure shows that oxygen-rich PNe 
possess the same spread in sulfur deficits as their carbon-rich counterparts. Furthermore, 
there is no correlation between the sulfur deficit and the occurrence of the "30 yU.m" feature; 
some carbon-rich PNe which exhibit the "30 yum" dust spectral feature also show low sulfur 
deficits or even enhanced sulfur abundances. Figures |9] and [10] show the spectral energy 
distributions for a selection of oxygen-rich and carbon-rich PNe, respectively. Although ISO 
SWS spectra are not available for all objects, it is clear that, for instance, NGC 7027 exhibits 
a "30 /im" feature but is barely sulfur- deficient (SD = 0.02); while IC 418, which is enriched 
in sulfur (SD = —0.25), also has a strong "30 /im" feature. 

If sulfur-bearing dust does not exhibit spectral features (only contributes to continuum) 
and can exist in an 0-rich environment, we would still expect to see a correlation between 
the dust abundance and the sulfur deficit. This is not the case. PNe typically show double- 
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humped SEDs, with a short-wavelength hump due to starhght and a long-wavelength hump 
due to dust emission. Examination of the SEDs presented in Figures [9] and [TO] show that there 
is no correlation between the S-deficit and whether starlight or dust emission dominates the 
SED. In order to quantify the ratio of starlight to dust emission we have determined the flux 
contained in the starlight and dust emission humps by interpolating between the photometry 
points and calculating the flux beneath each hump. The ratio of starlight to dust emission 
as a function of sulfur-deficit is plotted in Figure [TTl 

In some cases, the photometry measurements available are sparse (e.g. IC 2165), which 
makes the determination of the starlight to dust emission ratio uncertain. Consequently 
Figure [11] also contains the ratio of the V-band photometry to that in the infrared for the 
25, 60 and 100 /im IRAS photometry - as these far-IR photometry points should be good 
proxies for dust emission strength. We cannot find any correlation between dust emission and 
the sulfur-deficit and thus we rule out dust formation as the source of the sulfur abundance 
anomaly. 

Another potential sink for atomic/ionized sulfur gas is molecule formation. Like dust 
grains, molecules form during the AGB phase . Sulfur-bearing rnolecules observed aroun d 



AGB stars include SO, SO2, H2S, CS and SiS flOmont et al.lll993l : iBuiarrabal et al. Ill994l ). 
and can exist in both C-rich and 0-rich environments. As the star evolves towards the PN 
phase, the molecules will be destroyed by the increasingly high energy photons emanating 
from the star. Howeve r, molecules can survive in highly evolved PNe (e.g. the Helix Nebula; 



Matsuura et al.l [20091 . and references therein) as long as they are shielded from the high 
energy photons by dust (possibly in clumps). Sulfur-rich molecules will not return their 
sulfur bounty to the gas if the dust shields them. In this case we would still expect to 
see a correlation between the dust abundance and the sulfur depletion. As shown above 
(Figures [^ and [TUj) . this is not the case. However, dust emission in clumpy nebulae (like 
the Helix Nebula) is weak even though the presence of dust is inferred from the survival of 
molecular species. 

In order to test the hypothesis that sulfur can be sequestered in molecules in clumpy neb- 
ulae we sought correlations between the ratio of ionized/molecular mass in PNe with the sul- 



fur de fi cit ( Figure [T^. The r atios of molecular to ionized gas are taken from [Huggins fc Healy 



(Il989[ ): [Huggins et al.l ( ]1996[ ). Unfortunately, the number of PNe for which we have data on 
both molecular mass and sulfur anomolies is too small to be statistically significant. How- 
ever, the limited sample in Figure [12] suggests that the highest S-deficit is associated with 
the lowest molecule content. 



We have shown that sulfur abundance anomalies cannot be attributed to sequestration of 
sulfur into dust grains. Formation and protection of sulfur-bearing molecules may provide a 
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mechanism for sulfur sequestration. Further study of correlations between sulfur abundances 
and dumpiness of the host nebulae, and with evolutionary status of the nebulae are required. 



5. Nucleosynthesis in Asymptotic Giant Branch stars 



The illuminated PN is comprised of material from the deep convective envelope, hence 
nebular abundances should reveal information about the efficiency of mixing events and 
chemical processing that took place duri ng previous evolutionary phases, in addition to 



the initial compositi on of the parent star f lDopita et al. 1 119971 : iKarakas fc Lattanzio 1 12003 



Karakas et al. 1120091 ). In light of this information, is it possible to explain the sulfur anomaly 
in PNe by examining the observed or predicted nucleosynthesis outcomes from AGB stars? 
We start with a brief review of AGB evolution and nucleosynthesis. 



During the TP-AGB phase the He-burning shell bec omes thermally unst able e very 10 



years 



or so, depending on the core mass. We refer to IBusso et al. I (119991 ) and iHerwig 



( 120051 ) for reviews of AGB evolution and nucleosynthesis. The energy from the thermal 
pulse (TP) drives a convective pocket in the He-rich intershell, that mixes the products of He- 
nucleosynthesis within this region. The energy provided by the TP expands the whole star, 
pushing the H-shell out to cooler regions where it is almost extinguished, and subsequently 
allowing the convective envelope to move inwards (in mass) to regions previously mixed by the 
flash-driven convective pocket. This inward movement of the convective envelope is known 
as the third dredge-up (TDU), and is responsible for enriching the surface in ^^C and other 
products of He-burning, as well as heavy elements produced by the slow neutron capture 
process (the s process). Following the TDU, the star contracts and the H-shell is re-ignited, 
providing most of the surface luminosity for the next interpulse period. In intermediate- 
mass AGB stars with initial masses > 4M0 (with core masses > O.SMq), the base of the 
convective envelope can penetrate into the top of the H-burning shell, causing proton-capture 
nucleosynthesis to occur there (hot bottom burning, HBB). The surface composition of the 
star will be changed owing to the fact that the entire envelope is exposed to the hot burning 
region a few thousand times per interpu lse period. Given the fast evol utionary timescales 
of post-AGB stars with Mcorc ^ O.SM© (IBlocker &: Schoenberner Ill99ll ) it is not clear how 
many end up illuminated as PNe. For this reason, we will concentrate on nucleosynthesis in 
lower mass AGB stars. 

In low-mass AGB stars the evolution of the composition of the surface layers is primarily 
affected by the occurrence and efficiency of the TDU. The efficiency of the TDU and the 
composition of the He-intershell have been s hown to vary as a function of t he initial mass, 



metallicity, as well as time during the AGB (IBoothroyd fc Sackmann Ill988l : iKarakas et al 
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20021 : iStraniero et al. 1 120031 : iKarakas Il2010[ ). Observationally these quantities are not well 
known but a few general trends can be inferred. At the metallicity of the Galactic disk, the 
TDU likely occurs in stars with initial masses > I.5M0; this mass is reduced to ~ IMq at 
the metallicities of the Mage llanic Clouds [with metallicities ~ 2 — 5 times less than solar 
(jWallerstein &: Knapp Ill998l )]. Of particular interest for this study is the nucleosynthesis 
that occurs in the He-intershell. This is harder to determine observationally. The subclass 
of PG1159 post-AGB stars are H- deficient and show w hat has been considered to be He- 
intershell material at their surface (jWerner et al. II2009I ). In these stars, C abundances vary 
from 15-60% (by mass), and O from 2-20%. There are also trace quantities of other He- shell 
burning products including Ne (~2%, probably in the form of ^^Ne) and F. Phosphorous is 
solar, sulfur has abundances v arying from 0.01 tirnes so lar to solar, and some PG 1159 stars 
are also highly deficient in Fe (jWerner fc Herwig II2006I ). 



The C and O PG 1159 abundances are in direct contrast to intershell abundances 



predicted by stellar evolution models of AGB stars (e.g., iKarakas Il2010l ). Standard models 
predict that He constitutes about 75% (by mass), C roughly ~ 25%, and there is about 1—2% 
each of ^^Ne and ^^O. The p redicted abundances of F and Ne agree with the abundances 



measured in PG 1159 stars (jWerner &: Ranch 1 1 1994 IWerner et al. I l2005l ). Ne from the 



intershell is predicted to increase the surface Ne/0 ratio, but only when there is very efficient 
TDU and this has been pr e dicted to occur in on ly a narrow mass range (~ 2.5 — 3.5M0 



Karakas fc Lattanzio I l2003l : IKarakas et al. I l2009l ). The remarkably constant Ne/0 ratio 



observed in PNe with a wide variety of chemical compositions indicates that in most cases 
this ratio is not altered during the AGB. This is consistent with the conclusion that most 
PNe originate from masses < 2.5M0. Unlike C, O, and Ne, the abundances of P, S, and Fe 
are affected by neutron-capture reactions in the He-intershell. Depend ing on assumptions 



Busso et al. 2001 



made about t he efficiency of the ^^C(a,n)^^0 neutron source (see, e.g. 
Herwig II2005I ). P is predicted to be produced by factors of ~ 1.5 to 25 times the solar value 



whereas the intershell S abundances are 0.6 to 0.9 times solar depending on the AGB model 



flWerner fc Herwig II2OO6I : IKarakas et al. Il2009f ) 



To account for the discrepa ncy between t he observed PG 1159 abundances and the 
predicted intershell composition, iHerwig I ( 120001 ) adopted AGB models with overshoot into 
the C-0 core, that has a composition which is roughly 50% carbon and 50% oxygen. The 
inclusion of diffusive convective overshoot at the inner edge of the flash-driven convective 
pocket during a TP dredges up some of the C-0 ma terial and leads to large increases in the 
C and O abundances in the intershell (IHerwig II2OOOI ). in line with the PG 115 9 observations . 
The amount of overshoot applied is still however an uncertain free parameter (IHerwig II2OOOI ). 



Figs. [13] and [H] show the predicted abundances of oxygen, neon, and sulfur from the 
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AGB models [the Z 



from Karakas et al. 



0.02,0.008,0.004 models from iKarakas 
and the Z 



(2010): the Z = 0.01 model 



f l2010h : and the Z = 0.001 model is from lAlves-Brito et al. I fl201l[ )]. 
The surface abundances of the AGB models are taken at the tip of the AGB, after the 
last computed thermal pulse and are taken to be representative of the composition of the 
envelope when it is ejected from the star. In each figure we show 12+log(0/H) abundances 
on the X-axis versus 12+log(Ne/H) (Fig. [13]) or 12+log(S/H) (Fig. [14]). 

In Fig. [13] we see that there is considerable spread in the Ne abundances predicted from 
the AGB models. The models that produce the most Ne are t hose partial mixing zones an d 
those models with masses of 3Mq, which have efficient TDU fjKarakas fc Lattanzio II2003I ). 
Studies of the s-process in low-mass AGB stars suggest that protons need to be partially 
mixed into the top layers of the He-inte rshell in order to effic iently activate the main neutron 
producing reaction ^^C(a,n)^^0 (e.g., Gallino et al. IQQsi). In the AGB models shown in 
Fig. [13] and described for the most part in [Karakas I ( ]2010[ ). an exponentially decreasing 



proton profile is introduced into the top ~10-15% of the He-intershell (a partially mixed 
zone or PMZ). The protons are captured by the abundant ^^C to form a ^^C pocket (and 
usually a ^^N pocket), with neutrons released by the ^'^C(a,n)^^0 reaction. The extra Ne 
production arises from the extra ^^N, which is converted to ^^Ne during convective thermal 
pulses. From Fig[T3]we can see that the effect of this is small at higher metallicities or in lower 
masses (e.g., the I.8M0, Z = 0.01). The observational data for O and Ne (see Fig. [2]) shows 
some spread, and that most of the lower mass models (< 3M0) fall within that range. Initial 
mass function considerations, which favor lower mass stars, together with these predictions 
suggest that the majority Galactic disk PNe evolved from stars of < 2M0. There are a few 
PNe with higher Ne abundances and these may have evolved from higher mass progenitors 
of ~ 3M0. At the lowest metallicities shown in Fig. [13] we similarly conclude that most PNe 
with log(0/H)< 8 likely have evolved from progenitors less massive than I.5M0 and/or did 
not experience efficient third dredge- up. 

In Fig|Tl]we show the predicted abundances of O and S from the same set of AGB models 
shown in Fig. [13] In contrast to the predicted spread in Ne, there is no predicted spread in 
sulfur present in any of the models. We can conclude that AGB models do not result in a 
net production (or destruction) of S. This is consis tent with other AGB calculations (e.g., 
Forestini fc Charbonnel lllQQTt ICristallo et al. II2009I ) and suggests that current AGB models 
cannot solve the S anomaly problem. 

However, the lowest metallicity models show an increase of O and move to the right in 
the O versus S plot. What if all AGB models, regardless of metallicity, dredged- up sufficient 
O to move them to the right on the O versus S plot? This could be achieved if convective 
overshoot from the C-0 core increased the O content of the He-intershell as discussed above 
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for PG 1159 stars. In iKarakas et al. I (120101 ) the synthetic evolution of a 1.8Mq, Z = 0.01 
model with higher C and O intershell abundances was computed (it was assumed that the 
neon and sulfur abundance are not changed by this overshoot). In Figs. US] and [T3] we show 
the results of the synthetic I.SMq, Z = 0.01 AGB model with an O intershell abundance of 
20% (by mass), compared to the model with a standard intershell composition (the arrow 
connects the models and shows the direction caused by the shift in O). This model moves 
to the right of the line in both figures, and could provide an explanation for the S anomaly. 

Is there a need for such overshoot in AGB models? As indicated previously, plots 
of sulfur versus neon and sulfur versus argon for PNe compared to HII regions show an 
offset in the sulfur abundance by the same amount compared to the plot of sulfur versus 
oxygen. This indicates that the problem lies with the sulfur abundances and not with the 
AGB nucleosynthesis models. However, there are arguments for the inclusion of convective 
overshoot from the C-0 core into the envelope. These are largely based on the composition 
of the post-AGB PG 1159 s tars. These o bjects are He-rich, likely caused by a late or 
very late thermal pulse (e.g., [Blocker 1 1200 ll ). Another supporti ng piece of evidence is th e 
apparent correlation between oxygen isotope ratios in AGB stars (IHarris et al. lll985Ul987l ). 
Not only are they correlated with each ot her but appear to b e correlated with the observed 
carbon and s-process element abundances flHarris et al. 1119871 ). Given the large observational 
unc ertainties in the derived o xygen isotope ratios it is difficult to draw firm conclusions (see 
also lSmith fc Lambert Ill990l ). Lastly, the high ^^C intershell abundances that would result 
from overshoot into the C-0 core changes the neutron exposure and the resulting s-process 
abundance dist ribution such that th ey are no longer consistent with mainstream SiC grains 
measurements (ILugaro et al. II2003I ). This suggests that such overshoot is not common in 
AGB stars that become C and s-process rich. 

In summary, current AGB models are unable to account for the S anomaly observed in 
PNe. The inclusion of overshoot from the C-0 core into the He-intershell could provide a 
tentative explanation but only when considering the correlation between O and S. Overshoot 
does not help clarify the offset present between S and Ne and S and Ar. Furthermore, there 
are concerns as to the extent of convective overshoot into the C-0 during the AGB phase. 



6. Summary and Conclusions 

The sulfur anomaly is the tendency of planetary nebula sulfur abundances determined 
from coUisionally excited optical line measurements of and S"*"^ and an ionization correc- 
tion factor to fall significantly below expected levels. In order to understand its origin, we 
have quantified the sulfur anomaly by computing a sulfur deficit, i.e., the difference between 
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the observed and expected S abundances, and attempted to relate this value to a broad 
assortment of other factors relevant to PNe. These include directly measurable properties of 
the central star and nebula, problems related to determining the contribution of unobservable 
ionization stages of S to the total elemental abundance (including the effect of inappropri- 
ate rates of dielectronic recombination in nebular models used to determine the ionization 
correction factor for S), sequestration of S into dust and/or molecules, and nucleosynthesis 
in AGB stars which might either produce S-poor or 0-rich material. 

Nearly all of these factors can be eliminated from consideration as offering a solution to 
the problem. We find no evidence that the gap between observed and expected S abundances 
in PNe is statistically related to any of the 24 nebular or central star properties which we 
considered with the one exception of log(S"''^/S"'"^) which appears to exhibits a positive 
correlation. A careful comparison of sulfur deficits and possible IR signatures fails to reveal 
any indication of sequestration of sulfur into dust or molecules. Likewise for nucleosynthesis 
scenarios. Using AGB evolution models, we considered in detail if sulfur deficits could be 
related to either sulfur destruction or oxygen production during the AGB stage of stellar 
evolution. Here again, we obtained a null result. 

The remaining factors which we considered are associated with the sulfur ionization 
correction factor and the ability to correct for abundances of sulfur ions higher than S"^^. 
We presented photoionization models of three PNe representing a range in sulfur deficit. The 
models appeared to be consistent with true gas-phase abundance shortfalls, leaving us with 
the option of adopting that as the solution or speculating that the models (from which the 
ICFs are formulated) are failing to account for significant occupation of higher ionization 
stages. The latter seemed less radical, especially since our consideration of dust /molecules 
and AGB nucleosynthesis strongly argues against a truly reduced S abundance. Therefore, 
we further pursued the IGF angle. 

We evaluated the impact of under- or overestimated dielectronic recombination rates 
in altering the sulfur ion population distribution and thereby producing a spurious model- 
based IGF. We did this by computing large grids of photoionization models in which we 
tested factor-of-two changes in the rates and compared output line strength ratios with their 
observed counterparts. We concluded that no reasonable amount of uncertainty in DR rates 
could explain the sulfur deficits. 

Finally, we looked carefully at the PN sulfur abundances derived from direct observation 
of S^^ at 10.5/im. Having this measurement presumably reduces or eliminates the need to 
employ an IGF. We were unable to show with statistical certainty that the sulfur deficit for 
each object correlates with its S~^^ abundance. On the other hand, a visual comparison of 
positions in the e(S) versus e(0) plane of objects with S abundances inferred from the IGF 
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method and those whose abundances are derived with the IR method show that objects in 
the latter category generally have smaller sulfur deficits. In fact once the S~^^ abundances 
are accounted for, the average size of the remaining deficits for the IR objects resemble the 
expected size of the abundance scatter, although the deficits themselves are not completely 
eliminated. It would appear from this part of the analysis, then, that the sulfur deficit is 
related to the failure of ICFs commonly in use to correct for the true amount of S"*"^ (and 
possibly higher ionization stages) present in these objects. This notion is supported by the 
direct correlation we found in § 2 between log(S"'"'^/S"'"^) and the size of the sulfur deficit. 
Statistical confirmation of this should follow once additional IR observations are performed. 
Independent theoretical verification o f this point could come from model computations such 
as those in progress by IJacob I (120 111 ). By introducing a negative density gradient into his 
CLOUDY models, Jacob boosts the relative amounts of the higher ionization stages of sulfur. 
His approach stands in contrast to the usual use of a constant gas density over the entire 
nebula, i.e., the common assumption made when using models to infer the form of the sulfur 
ICF. 
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A. Correlation Coefficients 



When looking for correlations amongst various parameters it is important to define what 
constitutes a correlation. The strength of a correlation can be measured via the correlation 
coefficient (r). Often it is assumed that any linear regr ession fit for which r> 0.5 is an 



1997 



significant correlation as a set of points for which r > 0.7 (see also iDijkstra et al 



acceptable criterion for a significant correla t ion (s ee e. g.. ISpeck et al 
20051 ). At the same time [Thompson et al.l (120061 ) and iGuha Nivogi et al.l (1201 



Diikstra et al. 



) define a 



20051 ) 
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The ^sample correlation coefficient\ r, measures both the direction and strength of the 
relationship between an independent variable x and a dependent variable y. The value of the 
sample correlation coefficent is — 1 < r < 1 such that if r ~ there is no relationship between 
X and y. The sign (+ or — ) of r merely designates whether the relationship between x and 
y is positiv e or ne gative (i.e. it defines the direction of the slope of the trend.). According 



to 



Rumseyl (12003), most statisticians lik e to see |r| > 0.6 in order to consider a correlation 



significant, while 

a correlation. For the analysis here we choose the former, |r| > 0.6. 



Thompson et all (120061 ) stated r should be > 0.7, in order to be considered 
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Table 1. Sulfur Deficit Correlation Tests 
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^Variables in order from top to bottom are: reddening factor, log 
of the ion ratio, oxygen ionization correction factor, sulfur ionization 
correction factor, log of the nebular diameter, central star mass, log 
of the central star's luminosity, neon deficit (defined in the same way 
as the sulfur deficit), log oxygen abundance, log of the ion ratio (4 
rows), log of the difference between a PN's measured O abundance 
and the ambient interstellar abundance at its location along the disk, 
PN's galactocentric distance along the disk, electron temperatures 
derived from the lines indicated (5 rows), electron densities derived 
from the lines indicated (2 rows), and excitation class based upon the 
emission line indicated (2 rows). 



Table 2. Photoionization Models 
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Table 3. Comparison of ICF and IR Results 
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le(X) = 12 + log(X/H) 
2SVSoxp = 1 -log-i(-SDi,.) 

%haw et al. ] l|2010h . Values for e{Sicf) and SDi cf were derived from the and S^^ ion ic 
abundances listed in their paper and the sulfur ICF formulation of lKwitter fc Henry ] 1 I2OOII) . 

^Pottasch fc Bernard-Salas ] | |2010|'I and iHenrv et al. I 1 I2OO8I . DdDm-1 only). For all but 
DdDm-1, oxygen abundances in column 2 were taken from the HKB04 survey, while the sulfur 
abundances were those determined by Pottasch & Bernard-Salas. 
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Fig. 1. — Plot of e(S) versus e(0) for a large sample of MWG and M31 disk planetary nebulae 
(open symbols) along with a sample of H II regions and blue compact galaxies selected from 
the literature (H2BCG; filled symbols). The solid line is a least squares fi t to t he H2BCG 



sample. The dashed lines show solar abundances from lAsplund et al. I (120091 ). Average 



uncertainties for the S and O abundances of the PN sample, determined by propagating line 
measurement uncertainties through the abundance-determining step, are shown. 
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Fig. 2. — Same as Fig. [T] but for e(Ne) versus e(0). We have also added PNe from the LMC 
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Fig. 3. — Logs of the [S III]/[S II] and [O III]/[0 II] ratios plotted against each other for the 
same observed MW disk PN sample displayed in Fig. [T] (filled circles) and photoionization 
models using Atlas (open squares) or Ranch AGB (open circles) input SEDs. 



^ 1 5 

(N 

o 
+ 

(N . 
m 1 
m 

w 

O 



0.5 



• DRx 1 

o DRx 0.5 
o DRx 2 



yo o 




p.o 
,•• o* 



oo 



0.5 



1 1.5 

log[(5007+4959)/3727] 



2.5 



Fig. 4. — Logs of the [S ni]/[S II] and [O III]/[0 II] ratios plotted against each other 
for photoionization models only, using the standard DR rates (closed circles), one-half the 
standard rates (open diamonds), and twice the standard rates (open circles). 
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Fig. 5. — Sulfur abundances inferred from IR data versus S abundances derived using an 
ICF. Filled circles and diamonds show positions for PNe in the MWG disk and the SMC, 
respectively. The dashed line shows the one-to-one correspondence. 
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Fig. 6. — The difference in logarithmic S abundances inferred from IR and ICF methods 
versus SDjcf. Filled circles and diamonds represent PNe in the Milky Way disk and SMC, 
respectively. Objects to the right of the vertical dashed line have S abundances which fall 
below the values on the H2BCG line for their O abundance, while objects to the left of the 
dashed line have S abundances above that trend line. The horizontal dashed line separates 
the plot into regions of positive and negative values for Sir — Sicf- The solid line shows a 
least squares fit to the subset points with positive values for SDjcf- 




Fig. 7. — Same as Fi g . [T] but now showing t he objects f r om t he samples of 
Pottasch fc Bernard-Salas I (120 lOl . open diamonds) and IShaw et aL I (120 lOl . open squares) 
in which IR observations were used to infer S abundances. 
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(C/0) 

Fig. 8. — Sulfur deficit versus tlie observed C/0 abundance ratio for those sample objects 
for which carbon has been measured. The vertical dashed line divides the space between 
C/0<1 (oxygen-rich; triangles) and C/0>1 (carbon- rich; circles) environments. 
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Fig. 9. — Spectral energy distributions for oxygen-rich PNe. x axis is wavelength in /im; y 
axis is A Fa in Wm~^; filled triangles are the optical and near-IR photometry (collected from 
Simbad); filled squares are IRAS photometry; filled circles are MSX photometry; and solid 
line is the ISO SWS spectrum. Arrows indicate an upper limit on the photometry points. 
The S-deficit is indicated in the legend for each SED. 
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Fig. 10. — Spectral energy distributions for carbon-rich PNe. 

X axis is wavelength in /im; y axis is A Fa in W ; filled triangles are the optical and near- 
IR photometry (collected from Simbad); filled squares are IRAS photometry; filled circles 
are MSX photometry; and solid line is the ISO SWS spectrum. Arrows indicate an upper 
limit on the photometry points. The S-deficit is indicated in the legend for each SED. 
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Fig. 11. — The ratio of starlight to dust emission as a function of sulfur deficit. Left panel 
the ratio of starlight to dust emission is determined by calculating the area under each hump 
in the SED for each nebula. Right panel the ratio of starlight to dust emission is estimated 
using the V-band photometry as a proxy for starlight, and the IRAS 25, 60 and 100 /zm 
photometries as proxies for dust emission. In both cases x-axis is the ratio of starlight to 
dust emission; the dustiness of the nebulae increases to the left, y-axis is the S-deficit. 
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Fig. 12. — The ratio of the molecular gas mass to ionized gas mass as a function of sulfur 
deficit. 
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Fig. 13. — The predicted loge(0/H) versus loge(Ne/H) abundances from a selection of AGB 
models. The composition of the model is sampled from the surface of the star at the very 
tip of the AGB. The solar composition is shown, and the approximate relationship between 
O and Ne abundances from HIT regions (dotted line). Initial abundances used in the model 
calculations are shown by the cross symbols for each respective metallicity. The legend 
shows which symbols correspond to models of a given metallicity. Predictions include: 3M0, 
Z = 0.02 models with and without a partially mixed zone (solid black squares); a 1.8Mq, 
Z = 0.01 model with and without a partially mixed zone (solid green triangles); I.SMq and 
3Mq models of Z = 0.008 (blue diagonal crosses); I.5M0 and 3M0 models of Z = 0.004 (sohd 
aqua stars); and a 1.5Mq, Z = 0.001 model with and without a partially mixed zone (magenta 
circle with cross-hatch). In each case the upper symbol for the same mass corresponds to 
the model with a partially mixed zone, as these have the highest Ne abundances. We also 
include results of the synthetic AGB model with an O intershell abundance of 20% (by mass), 
connected by the arrow to the model with a standard intershell composition (0.8% by mass), 
see the text for details. 
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Fig. 14. — Same as Fig. [13] but for the predicted loge(0/H) versus loge(S/H). Here there 
is no clear distinction between models of different mass or for models with and without a 
partially mixed zone, except at Z = 0.001 (where we only include results at I.5M0). Also 
included are results of the synthetic AGB model with an O intershell abundance of 20% (by 
mass), connected by the arrow to the model with a standard intershell composition (0.8% 
by mass); see text for details. 



